Activated sampling in complex materials at finite temperature: the properly-obeying-probability activation-relaxation technique

While the dynamics of many complex systems is dominated by activated events, there are very few simulation methods that take advantage of this fact. Most of these procedures are restricted to relatively simple systems or, as with the activation-relaxation technique (ART), sample the conformation space efficiently at the cost of a correct thermodynamical description. We present here an extension of ART, the properly-obeying-probability ART (POP-ART), that obeys detailed balance and samples correctly the thermodynamic ensemble. Testing POP-ART on two model systems, a vacancy and an interstitial in crystalline silicon, we show that this method recovers the proper thermodynamical weights associated with the various accessible states and is significantly faster than MD in the diffusion of a vacancy below 700 K.Comment: 10 pages, 3 figure

Algorithms for 3D rigidity analysis and a first order percolation transition

A fast computer algorithm, the pebble game, has been used successfully to study rigidity percolation on 2D elastic networks, as well as on a special class of 3D networks, the bond-bending networks. Application of the pebble game approach to general 3D networks has been hindered by the fact that the underlying mathematical theory is, strictly speaking, invalid in this case. We construct an approximate pebble game algorithm for general 3D networks, as well as a slower but exact algorithm, the relaxation algorithm, that we use for testing the new pebble game. Based on the results of these tests and additional considerations, we argue that in the particular case of randomly diluted central-force networks on BCC and FCC lattices, the pebble game is essentially exact. Using the pebble game, we observe an extremely sharp jump in the largest rigid cluster size in bond-diluted central-force networks in 3D, with the percolating cluster appearing and taking up most of the network after a single bond addition. This strongly suggests a first order rigidity percolation transition, which is in contrast to the second order transitions found previously for the 2D central-force and 3D bond-bending networks. While a first order rigidity transition has been observed for Bethe lattices and networks with ``chemical order'', this is the first time it has been seen for a regular randomly diluted network. In the case of site dilution, the transition is also first order for BCC, but results for FCC suggest a second order transition. Even in bond-diluted lattices, while the transition appears massively first order in the order parameter (the percolating cluster size), it is continuous in the elastic moduli. This, and the apparent non-universality, make this phase transition highly unusual.Comment: 28 pages, 19 figure

Self-organization with equilibration: a model for the intermediate phase in rigidity percolation

Recent experimental results for covalent glasses suggest the existence of an intermediate phase attributed to the self-organization of the glass network resulting from the tendency to minimize its internal stress. However, the exact nature of this experimentally measured phase remains unclear. We modify a previously proposed model of self-organization by generating a uniform sampling of stress-free networks. In our model, studied on a diluted triangular lattice, an unusual intermediate phase appears, in which both rigid and floppy networks have a chance to occur, a result also observed in a related model on a Bethe lattice by Barre et al. [Phys. Rev. Lett. 94, 208701 (2005)]. Our results for the bond-configurational entropy of self-organized networks, which turns out to be only about 2% lower than that of random networks, suggest that a self-organized intermediate phase could be common in systems near the rigidity percolation threshold.Comment: 9 pages, 6 figure

Self-organized criticality in the intermediate phase of rigidity percolation

Experimental results for covalent glasses have highlighted the existence of a new self-organized phase due to the tendency of glass networks to minimize internal stress. Recently, we have shown that an equilibrated self-organized two-dimensional lattice-based model also possesses an intermediate phase in which a percolating rigid cluster exists with a probability between zero and one, depending on the average coordination of the network. In this paper, we study the properties of this intermediate phase in more detail. We find that microscopic perturbations, such as the addition or removal of a single bond, can affect the rigidity of macroscopic regions of the network, in particular, creating or destroying percolation. This, together with a power-law distribution of rigid cluster sizes, suggests that the system is maintained in a critical state on the rigid/floppy boundary throughout the intermediate phase, a behavior similar to self-organized criticality, but, remarkably, in a thermodynamically equilibrated state. The distinction between percolating and non-percolating networks appears physically meaningless, even though the percolating cluster, when it exists, takes up a finite fraction of the network. We point out both similarities and differences between the intermediate phase and the critical point of ordinary percolation models without self-organization. Our results are consistent with an interpretation of recent experiments on the pressure dependence of Raman frequencies in chalcogenide glasses in terms of network homogeneity.Comment: 20 pages, 18 figure

Application of microfluidic systems in modelling impacts of environmental structure on stress-sensing by individual microbial cells

Environmental structure describes physical structure that can determine heterogenous spatial distribution of biotic and abiotic (nutrients, stressors etc.) components of a microorganism's microenvironment. This study investigated the impact of micrometre-scale structure on microbial stress sensing, using yeast cells exposed to copper in microfluidic devices comprising either complex soil-like architectures or simplified environmental structures. In the soil micromodels, the responses of individual cells to inflowing medium supplemented with high copper (using cells expressing a copper-responsive pCUP1-reporter fusion) could be described neither by spatial metrics developed to quantify proximity to environmental structures and surrounding space, nor by computational modelling of fluid flow in the systems. In contrast, the proximities of cells to structures did correlate with their responses to elevated copper in microfluidic chambers that contained simplified environmental structure. Here, cells within more open spaces showed the stronger responses to the copper-supplemented inflow. These insights highlight not only the importance of structure for microbial responses to their chemical environment, but also how predictive modelling of these interactions can depend on complexity of the system, even when deploying controlled laboratory conditions and microfluidics

Droplet coalescence is initiated by thermal motion

The classical notion of the coalescence of two droplets of the same radius R is that surface tension drives an initially singular flow. In this Letter we show, using molecular dynamics simulations of coalescing water nanodroplets, that after single or multiple bridges form due to the presence of thermal capillary waves, the bridge growth commences in a thermal regime. Here, the bridges expand linearly in time much faster than the viscous-capillary speed due to collective molecular jumps near the bridge fronts. Transition to the classical hydrodynamic regime only occurs once the bridge radius exceeds a thermal length scale lT∼√R